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F. INTRODUCTION 


At the present time the sphere represents the only radar target 
of finite extent for which numerical scattering results may be obtajned 
directly and simply from the rigorous eigenfunction solution of the plane 
wave scattering problem. The rigorous eigenfunction solution of the 
thin metallic disk problem has been available in the literature for a 
considerable period of time [1]; however, only limited numerical results 
have become available due to the difficulty of the numerical compu- 
tations required. 


This state of affairs is quite unfortunate since the disk offers 
significant advantages as a standard radar calibration target when com- 
pared with the sphere. First, precision machining of a disk is much 
simpler than that of a sphere; and, second, precise alignment of the 
target for phase measurements is considerably simpler for a disk than 
for a sphere. 


Furthermore, a great deal of basic understanding of scattering 
mechanisms is potentially available from the study of the thin disk. 
This is a consequence of the fact that it is the only target of finite 
extent, other than the sphere, for which a rigorous solution is avail- 
able; and it is the only case for targets having a sharp edge. Thus, 
a complete understanding of scattering by a disk would provide another 
canonical solution to complement that of the sphere. 


For these reasons, earlier work at The Ohio State University Electro- 
Science Laboratory [2,3,4] has been extended to generate a computer program 
capable of handling the general problem of far field scattering of a 
plane wave by a thin, metallic disk. This program permits incident plane 
waves of arbitrary incidence direction and arbitrary polarization and 
computes the amplitude and phase of both components of the scattered 
field at any point on the far field sphere. 
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The program has been generated in a user oriented form so that 
it can readily be used by any investigator without a detailed knowledge 
of the program. The program requires about 16K of core memory and executes 
in a matter of seconds on the ESL Datacraft 6024 computer operating in 
a time sharing mode. 


II. THE GEOMETRY 


The disk of radius a is centered at the or“gin and lies in the 
x-y plane. The direction of propagation of the incident plane wave is 
taken to be in the x-z plane without loss of generality. The angle of 
incidence, Q» js measured from the positive z-axis, i.e., the normal 
to the disk, as shown in Figure 1. The scattered far field is to be 
evaluated in a direction specified by the conventional spherical coordin- 
ates, 8, and o.- 


The polarization of the incident E-field is aligned at an angle 
of a measured from the plane of incidence in the + direction. Thus, 
a=0 is associated with the parallel, E-plane, or 9-polarized case, and 
a= 3 is associated with the perpendicular, H-plane, or $-polarized case. 
Both the 6 and @ components of the scattered E-field will be determined. 


III. THE EIGENFUNCTION SOLUTION 


The solution presented here parallels that obtained by Andrejewski 
[1]. For convenience in the computation, the solution has been cast 
in terms of trigonometric rather than exponential functions. And, the 
more conventional notation of Flammer 15] has been followed. 


The incident electric field intensity is given by 


i a 
E ES (-cos@,cosaa, 


+ a 
sinaa,, 


i- 
+sin® cosa, )e!(K -Ptut) (1) 


2 


) 
' 
} 
} 


Figure 1. The geometry of the scattering problem. 
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Then, the scattered electric field intensity at an arbitrary point, 
(750.566) on the far field sphere may be expressed as 


iE 
ES = 552 ($254. G+ sina &,) ekrmat), (5) 
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It is convenient to express the scattered field E® in terms of a normalized 
field ES as 


E : 
ES = = © e7 i(kr-wt) - (6) 


This choice of normalization yields a particularly simple form for the 
normalized radar cross section, i.e., 
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The normalized electric field intensity in this case is given by 


(8) 
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The functions Y.. are given in terms of the spheroidal radial functions, 
R(')(4¢340), and the spheroidal angular functions, Smn(-ic,cosé), by 
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The prime on the summation symbol emphasizes the fact that the summation 
over n proceeds by increments of 2 as a consequence of the condition 
that n-m is even. The normalization function, Nan? and the spheroidal 
functions will be described later. 


The U and X functions are given by 
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Finally, the W and wp functions are given in terms of the spheroidal 
functions by 
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The angular spher::'dal functions may be expressed in terms of con- 
ventional spherical Legendre polynomials, prer(COS6) » as [5] 


Smn(7icscose) = % ann (-ic)PM, (cose) (16) 


where the prime indicates that the summation is over even values of r 

if n-m is even and over odd values of r if n-m is odd. The expansion 
coefficients, dmn(-ic), are those used by Flammer [5] and may be computed 
readily using a technique described in Reference 2. 


Since only spheroidal radial functions of zero argument and n-m 
even are required, the special relationships for these cases as presented 
by Flammer may be used: 
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The expansion coefficients, ann(-ic), used here are identical to those 


used in Equation (16). The radial spheroidal functions of the qth kind 
are found simply as in the spherical case by 


R19) (L4c340) = RO) (~4es40) - an?) (Ges 40). (19) 
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The normalization functions, Nan’ ic) are those required to cause 
the angular spheroidal functions to satisfy 


Santis) = PH(o) (20) 


and are given by : 
©, (r+2m)!Ld"(-ic)] 
Nan (7 1c) =2 toy —Tartamttyri ° (21) 


Equations (8) through (21) provide the complete solution to the 
general far field scattering problem. 


IV. THE PROGRAM 


Using the solution given in the preceeding section, a Fortran com- 
puter program was prepared for the calculation of far field scattering 
by a circular metallic disk. The program was developed on a time-sharing 
Datacraft 6024 machine having a 24 bit word length. The storage require- 
ments are 


Main Program and Subroutines (4,790) 15 words 
Library Subroutines (4,514) 15 words 
Common Storage (7,300) 15 words 
Total Storage (16,604), words 


(not including operating system and 1/0 buffers). 


This program executes a complete far field scattering computation in 
a matter of seconds in the time-sharing environment. 


A sample of the program I/0 is shown in Figure 2. In general, 
one need only provide 8 variables for the initial case to be executed; 
they are: 


1. KA = the electrical radius of the disk 
THETA INCIDENT = 6, Ldegrees] (see Figure 1) 
POLARIZATION = a [degrees] 

- THETA SCATTERED = 9, [degrees] 

- PHI SCATTERED = 9, [degrees] 
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WHICH VARIABLE IS TO BE INCREMENTED? 
if 1 thru 5: then the variable associated with 
that index above will be incremented 
if not 1 thru 5: then only the initial case will be 
calculated. In this event the remaining 
parameters are not requested. 


TYPE NUMBER OF CASES; enter the total number of computations 
to be performed. 


WHAT IS THE INCREMENT: enter the increment by which the variable 
selected should be increased after each execution. 


All inputs may be made in free format. 


The program labels the is column with the name of the variable 
to be incremented. The gre and 3rd columns contain the cross sections 
(Equation (7)) associated with the 9 and » components of the scattered 
field. The last four columns list the magnitudes and phases (in degrees) 
of both the @ and » components of the normalized scattered electric field, 
5 (Equation (8)). It should be noted that all of the elements of the 
scattering matrix are available in lines 95-104 of the program. 
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(TYPE “ESC™ TO RESTANT PPCSAM? 

(TYPE KABA TO STOP PROGRAM) 
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Figure 2. I/0 listing. 
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Upon completion of the cases desired, the program will request 
a new disk size and proceed as before. At any time one may type "ESC" 
and cause the current task to be terminated; the program will then again 
request a new disk size and proceed as before. 


Entering a disk size of 0 will cause the program to terminate. 
Entering a disk size of -1 will cause a brief statement of the problem 
geometry to be printed; following this a new disk size will be requested. 


A simplified flow diagram showing the computational logic is presented 
in Figure 3. The logic is quite straight forward and proceeds along 
the line specified by Equations (7-21). 


It is necessry to compute the eigenvalues of the spheroidal wave 
equation, Amn(-ic), in order to determine the expansion coefficients 
dr(-ic)s appearing in Equations (16), (17), (18) and (20). The eigenvalues 
are computed by the bisection method and the expansion coefficients are 
computed by recursion as described in Reference 2. 


The solution of the scattering problem consists of a triple summation 
over the indices m, n, and r. The truncation of these summations is 
performed internally by the software. Various functions are examined 
to determine if they are near the machine overflow level, i.e., 10°8 
for the Datacraft 6024. If this level is reached, the appropriate sum- 
mation is truncated as described in Appendix B. This procedure yields 
the best possible convergence for a machine having this dynamic range. 

This procedure has also been successfully used for sphere scattering 
calculations. In both cases the truncation is controlled largely by 
the tendency of the radial functions appearing in denominators to become 


extremely large. This tends to insure adequate convergence of the 
solution. 
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Figure 3. 
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Simplified flow diagram. (The numbers in parentheses 
refer to line numbers in the program listing presented 
in Appendix I). 
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SAMPLE RESULTS 


Some sample results are included here to serve as check cases and 
to demonstrate the utility of the program. 


Calculations were made as a function of electrical disk size for 
normal incidence backscatter. These results are tabulated in Figure 2; 
and the normalized radar cross section and scattered E-field phase are 
plotted in Figures 4 and 5, respectively. The Rayleigh or low frequency 
approximation 


cS, B(ka)® 
MRay 1 
and the Geometrical Theory of Diffraction or high frequency approximation 
, » (tka + 3H), > -i(4ka - §) 
— = — e - Ika * Geka ° - ika 


are also shown presented in these figures for comparison. 


Results for normal incidence, bistatic scattering are shown in 
Figures 6 and 7 as a function of the scattering direction. Both E- and 


H-plane results are given here for a disk size of ka=10 (diameter = 3.18). 


Results for specular bistatic E- and H-plane scattering from a 
disk of ka=10 are shown in Figures 8 and 9. In this case 6,70, and 
@, =180°. Note that the phases for 6.=0 in Figures 7 and 9 differ 
by 180°; this is a consequence of the chosen coordinate system. The 
calculations in Figure 7 were done for ,=0 and those in Figure 9 for 
$,=180°. 
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Figure 4. 
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Normal incidence, backscatter. 
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Figure 5. Normal incidence, backscatter. 
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Figure 6. Normal incidence, bistatic scatter. (ka=10). 
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Figure 7. Normal incidence, bistatic scatter. (ka=10). 
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Figure 8. Bistatic specular scatter. (ka=10) 
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Figure 9. Bistatic specular scatter. (ka=10) 


VI. | SUMMARY 


The expressions for calculating far field scattering characteristics 
of a plane wave by a thin metallic circular disk were presented. These 
expressions are applicable for any incident polarization or direction 
of incidence and for both scattered polarizations at any point on the 
far field sphere. This formulation was used to develop a user oriented 
computer program which is also described herein. This program has been 
used successfully for disk sizes up to ka=15. It yields both the radar 
cross section and phase of both scattered field components. The program 
executes in a matter of seconds and requires about 16k of core memory. 
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APPENDIX I 
PROGRAM LISTING 


The following is a Fortran listing of the program described in 
the body of this report. The use of this program is described in the 
section entitled the Program. Comments on the program are included in 
Appendix II. 
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APPENDIX ITI 
PROGRAM NOTES 


The following notes are intended to clarify the program listed 
in Appendix I. The notation LN will be used to denote the line number 
in the listing, i.e., the first number on each line. 


The first 9 lines of this program are non-executable and simply 
establish variable types and dimensions. LN1IO through 69 accept input 
data and initialize variables. LN2, 20, 150 and 151 are associated with 
the interactive 'ESC' which permits interruption of the program and return 
to the start of the data input segment. 


As noted in the body of the report, the scattered field is expressed 
as a triple summation over indices m, n, and r. In general, these indices 
range from 0, m, or 1 to some maximum value. In this program these indices 
are usually replaced by MM, NN, and IRR, respectively, where 


MM = m+] 
NN = (n-m)/2+1 
IRR = 1,2,3,... ‘ 


Thus, for computational convenience, these indices all range upward from 
1 in integer steps. 


The truncations of these summations are determined on the basis 
of tests of key variables. During execution these variables are tested 
and upon reaching absolute values in the range 1030 to 10°8, the appropriate 
summation is truncated. The line numbers associated with these teSts 


and the subsequent actions are tabulated below: 


* 
Aen eh nse AS TE NR bt Ye pmb 


m n r 


LN72 LN236 LN290 

78 257-259 294 

91 296-303 305-306 
388 331-332 327-330 
410-411 360-361 

434-436 364-365 

539-541 413 

551-552 


The necessary functions are formed by subroutine calls in Loop 
10 and Loop 13 (LN70-94). These functions are then combined to form 
the scattered E-fields and cross sections in LN95-128. LN130-132 provide 
the program output; and LN134-136 provide incrementing of the variable 
desired. LN138-149 contains the description of the coordinate system 
to be provided for ax input of KA=-1. The remainder of the listing con- 
sists of the required function subroutines. 


Whenever possible, variable names associated with the symbols presented 
in this report are used. One exception is the introduction of 


gorse ec) er 
Nan(~C)Rin (-ic;io) 


Other exceptions are the Z functions appearing in LN95-98 and LN554-595. 
These are functions of U, X and Y found in Equation (9) of the text. 


